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Abstract. We present an analytic description of the r- 
mode instability in newly-born neutron stars, using the 
approximation of uniform density. Our computation is 
consistently accurate to second order in the angular ve- 
locity of the star. We obtain formulae for the growth-time 
of the instability due to gravitational-wave emission, for 
both current and mass multipole radiation and for the 
damping timescale, due to viscosity. The I = m = 2 
current-multipole radiation dominates the timescale of the 
instability. We estimate the deviation of the second order 
accurate results from the lowest order approximation and 
show that the uncertainty in the equation of state has only 
a small effect on the onset of the r-mode instability. The 
viscosity coefficients and the cooling process in newly-born 
neutron stars are, at present, uncertain and our analytic 
formaulae enables a quick check of such effects on the de- 
velopment of the instability. 
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1. Introduction 

The recently discovered r-mode instability 
(Andersson 1998) in rotating neutron stars, has sig- 
nificant implications on the rotational evolution of a 
newly-born neutron star. The r-modes are unstable due 
to the Chandrasekhar-Friedman-Schutz (CFS) mecha- 



nism dChandrasekhar 1970| , Friedman & Schutz 1978) (see 
also Friedman & Morsink 1997). Andersson, Kokkotas & 
Schutz (1998) (see also Lindblom et al. 1998) find that 
the r-mode instability is responsible for slowing down 
a rapidly rotating, newly-born neutron star to rotation 
rates comparable to that of the initial period of the Crab 
pulsar (~19 ms) or the recently discovered 16 ms X-ray 
pulsar in the supernova remnant N157B (Marshall et 
al. 1998) (with an estimated initial period of 6-9 ms). 
This is achieved by the emission of current-quadrupole 
gravitational waves, which reduce the angular momentum 



of the star. Additionally, as the initially rapidly rotating 
star spins down, an energy equivalent to roughly 1% of 
a solar mass is radiated in gravitational waves, which 
makes the process an interesting source of detectable 
gravitational waves (Owen et al. 1998). 

In the present paper, we investigate the r-mode insta- 
bility to 2nd order accuracy in the angular velocity of the 
star, in the approximation of uniform density (the actual 
density profile of realistic neutron stars is nearly uniform) , 
where all results can be obtained analytically. |^ While our 
analytic results provide an independent check of the nu- 
merical results in Andersson et al. 1998 and Lindblom et 
al. 1998, our main objective is to present a simple set of 
equations, which enable one to obtain a qualitative in- 
sight into the mechanism of the r-mode instability and to 
quickly check the dependence of the instability on various 
important factors, such as the central density of the star, 
the different types of viscosity in neutron stars, the differ- 
ent possible cooling processes etc. Additionally, we expect 
that there is a number of issues related to the spinning 
mechanisms of pulsar such as accretion disc induced spin 
up, or the creation of millisecond pulsars due to accretion - 
induced collapse of a white dwarf ( Andersson et al. 1998 ), 
for which one can use the simple relations provided here 
for a fast but still accurate evaluation of the various evo- 
lution scenarii. 



2. The r-mode Instability 

Oscillations of stars are commonly described by the La- 
grangian displacement vector ^, which describes the dis- 
placement of a given fluid element due to the oscillation. 
Since ^ is a vector on the {9, (p) 2-sphere, it can be ana- 
lyzed into a sum of spheroidal and toroidal components (or 
polar and axial components, in a different terminology). 



The computation in Andersson et al. 1998 



is carried out in 
the 2nd order accurate, s low-ro tation formalism for N ewtonian 
stars by Saio ( Saio 1982 ) while Lindblom et al. 1998 , arrive at 



the same conclusion, using a first-order, approximate, calcula- 
tion. 
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In a non-rotating star, the usual /, p and g modes of oscil- 
lation are purely spheroidal, characterized by the indeces 
(Z, m) of the spherical harmonic function . In a rotating 
star, modes that reduce to purely spheroidal modes in the 
non-rotating limit, also acquire toroidal components. Con- 
versely, r-modes in a non-rotating star are purely toroidal 
modes with vanishing frequency. In a rotating star, the 
displacement vector acquires spheroidal components and 
the frequency in the rotating frame, to first order in the 
rotational frequency 17 of the star, becomes 



UJr 



2mn 
l{l + l) 



(1) 



for a given (l,m) mode. An inertial observer, measures a 
frequency of 



mf2 



(2) 



From (H) and (||), it can be deduced that, a counter- 
rotating (with respect to the star, as defined in the co- 
rotating frame) r-mode, appears as co-rotating with the 
star to a distant inertial observer. Thus, to 0(57), all r- 
modes with I > 2 are generically unstable to the emis- 
sion of gravitational radiation, due to the Chandrasekhar- 
Friedman-Schutz (CFS) mechanism (note that the I = 1 
r-mode is marginally unstable, to this order). The insta- 
bility is active for as long as its growth-time is shorter 
than the damping-time due to the viscosity of neutron- 
star matter. Its effect is to slow-down, within a year, a 
rapidly rotating neutron star to slow rotation rates and 
this explains why only slowly-rot ating pulsars are asso - 
ciated with supernova remnants ( Andersson et al. 1998|) . 
Thus, the r-mode instability does not allow millisecond 
pulsars to be formed after an accretion-induced collapse 
of a white dwarf. It seems that millisecond pulsars can 
only be formed by the accretion-induced spin-up of old, 
cold, neutron stars. 



3. The 2nd-order Accurate Slow Rotation 
Formalism 

To 0(il), the star is still spherical, and one can only de- 
termine the angular dependence of the r-mode eigenfunc- 
tions and their lowest order r-dependence (the latter is 
obtained by taking the curl of the perturbed equations 
of motion). For obtaining the second-order correction to 
the eigenfunctions and to the frequency, one must pro- 
ceed to a consistent 0(i7^) calculation. We follow the 
formali sm for com puting r-modes in Newtonian stars, 
that was presented in more detail in 



due to 3aio 1982 



Andersson et al. 199S . Here we will only summarize the 
equations needed for the uniform density case. 

3.1. Assumptions 

We make the following assumptions: 



1. the perturbations are adiabatic, 

2. the star has uniform density, 

3. the rotation of the star is uniform, and 

4. the perturbation of the gravitational potential can be 
neglected (Cowling approximation). 

These assumptions are justified by the fact that, even for 
temperatures T — lO^/iT, the thermal energy of the star 
is much less than the Fermi energy of its interior (> 60 
MeV). Also, at such temperatures, the initially differen- 
tially rotating proto-neutron star is ro tating uniformly . 



due to the formation of a solid crust (see Stergioulas 1998 



for a recent review on rotating neutron stars). The Cowl- 
ing approximation has been shown to yield sufficiently ac- 
curate results for r-modes in slowly rotating, Newtonian 



stars (3aio 1982, Provost et al. 1981) 



3.2. Definitions 

In a slowly rotating star, the dominant correction to its 
structure is of O(ri^). The analysis of perturbations of the 
star is simplified by introducing a new radial coordinate 
a, defined through 

r = a(l + e) , (3) 

where e = e(a, 0) is a quantity of 0(f2^), representing the 
deformation of the equilibrium structure from the non- 
rotating configuration. In the new coordinate system, all 
equilibrium quantities are functions of a only and the sur- 
faces of constant a are equipotential surfaces. 

Since equilibrium neutron stars are stationary and ax- 
isymmetric, a general oscillation can be analyzed into a 
sum of normal modes, with harmonic time-dependence 
gi(m0+wit) rpj^g displacement vector for a given r-mode can 



be written as (Saio 1982) 
i/a = T + S 



E 



H^^j_ sin" 



'ed^)Y>f, (4) 



where T and S are the toroidal and spheroidal parts of 
the displacement, respectively. Note that the toroidal part 
has vanishing a-component and is described only by the 
function Kim = Kim{o,), which multiplies a toroidal an- 
gular vector. The spheroidal part has an non- vanishing 
a-component of 0(f2^), described by the functions S^^ — 
Svii {o) ■ The and 0-components of the spheroidal part are 
described by the functions Hi,^^ — Hy^{a) (also of 0(fi^)), 
multiplying spheroidal angular vectors. 

The perturbation in the pressure is expressed in terms 
of spheroidal radial functions Ci^^t, as 

5p^ pgaY,(:,^Yl\ (5) 

where 5p is the Eulerian variation in the pressure (the 
variation of the pressure at a fixed point in space), p is 
the density and g = —p^^dP/da is the acceleration of 
gravity. 
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3.3. The Propensity Rule 

For zero-temperature (barotropic) stars, it can easily be 
shown from the perturbation equations, that only modes 
with I — ra exist. Then, only the v = I + \ , ji ^ m terms 
contribute in the expansions for the displacement vector 
and the perturbation in the pressure (in the remainder of 
the text, we will drop the index m in these quantities). 
The absence oil — 1 terms (the spherical harmonics 
are zero) means that rotation excites only higher multipole 
spheroidal parts. This is in agreement with the "propen- 
sity" rule suggested by Chandrasekhar & Ferrari 1991, for 



the oscillations of slowly rotating relativistic stars, i.e. 
that the rotational coupling of a toroidal Z-term with a 
spheroidal I -\- 1-terms is strongly favored over the coupling 
with a spheroidal I — 1-term. We find that, in uniform den- 
sity stars, the "propensity" rule completely eliminates the 
coupling to lower-multipole terms. 



In (0), h is 



(2/ + 3)— +6(/- l)i)2 

CTO 



+3 ?,D2 + a- 



da 



(14) 



Note that, the perturbation in the pressure is indepen- 
dent of the displacement vector and can be found by an- 
alytically integrating (|l0|), while the toroidal function Ki 
is given algebraically in terms of the perturbation in the 
pressure. The spheroidal function S';+i satisfies a differ- 
ential equation that depends on the perturbation in the 
pressure and the structure of the star and can not be ob- 
tained analytically, but will not be needed for the remain- 
der of this paper. The spheroidal function Hi is given al- 
gebraically in terms of Si-^-i and O+i- 



3.4. The Perturbation Equations 

A normal-mode solution to the perturbation equations 
satisfies the perturbed Euler equations, the perturbed con- 
tinuity equation and the relation between the perturba- 
tions in density and pressure. We define dimensionless fre- 
quencies as 

<gm) ' 



and 



\l/2 



f— ) 

\gmJ 



and expand the frequency in the rotating frame as 
UJr — (Ton? + a2'^'^ . (8) 
Writing the distortion parameter e as 

.2 



e = 



Di{a) + D2ia)P2icose) 



(9) 



(where P2 (cos 6) is the Legendre polynomial) and expand- 
ing the perturbations equations consistently to second or- 
der in the angular velocity of the star, we find that the 
eigenfunctions Q+i, Si+i, and Ki are given by the 

following set of equations: 



a^ = (/-l)Cm, 
da 

Hi+i = Si+i + 
and 



(4 + I) Si+i - hCi+i 
(2/ + l)(/ + l) 



(10) 
(11) 



8lao 



Ki 



.(/ + 1)^/21+3 
^ 



{l + l)a2+QD2 0+1,(12) 



(13) 



3.5. Boundary Conditions 

From the leading terms of Si and Q near a = 0, one obtains 
the boundary condition at the center of the star: 



{21 + 3)5i+i + hCi+i = 



(15) 



At the surface of the star, the Lagrangian variation of the 
pressure vanishes (a fluid element on the surface of the 
unperturbed configuration must also be on the surface of 
the perturbed configuration): 



(7) ^p = 6p + iVp = 0, 



or 



0+1 — Si 



i+i- 



(16) 



(17) 



To 0{n), (|T|) is satisfied trivially, while to 0{Vt^) it yields 
the correction to the eigenfrequency to that order. 



4. Eigenfunctions and Eigenfrequencies 

Eq. ( |l0|) for 0+1 implies a solution of the form 

0+1 - (18) 

Since O+i is of order 0{w'^), we normalize it to the di- 
mensionless quantity 



i-i 



^'+^ = "^ [r. 

Then, (pi) yields 



.(? + 1)^21+3 / f7w a 
^' = * 21 [-J[r) ' 



(19) 



(20) 



where R is the radius of the star. These are the only two 
eigenfunctions needed for the remainder of the paper. 
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An expression for the eigenfrequency of a given mode 5.2.1. Mass Multipoles 

can be obtained from constructing an integral relation us- mi i • i • i • r 7^ t • 

^, i 1 1 1 / ii T r ilie dommant mass multipole moment is oDi+i. in uni 

mg the perturbed iLuler equations (see the appendix of . \. r ■ ■ ■ c\.a- 



Saio 1982). For uniform density stars 
2 . 5/(^ + 1)2- 



o-Q = 
or 



l + V 
2 



and 



CT2 



{i + iy 



5?(/ + l)2-8 



(21) 
(22) 



form density stars, the Lagrangian variation of the density 
vanishes, Ap =0. From the relation between Lagrangian 
and Eulerian perturbations of a scalar quantity, it follows 
that 

5p = -£,Vp. (29) 
The derivative of the density across the surface is a Dirac 



l + l il + l)^ 

I c J- TT fRT : — : — i — 1 noi l\ T4- u ij u i J delta-function at a = i?, thus 

(ct. appendix 11 of Provost et al. 1981). It should be noted ' 

that the uniform density solution in appendix II of Provost 
et al. 1981 contradicts an erroneous conclusion in ap- 
pendix I. of the same paper. The authors conclude that The mass-multipole moment becomes 
there are no r- modes to OiVl^) in a barotropic star, by 
showing that the perturbation equations do not admit so- 
lutions with a Y[^{cos9) dependence of ^r- But the lowest- 
order angular dependence of is Yi^i (cos 9) . 



6p = ~aSi+iYi_^^pS{a - R) 
nultipc 
m^pR' 



i+i 



(30) 
(31) 



and, being of 0{i^^), it contributes to dE/dt\g 
0(f72'+8) term. 



5. Dissipation Time-scales 
5.1. Energy of mode 

The energy of the mode, measured in the rotating frame, 



IS 

E 



5.2.2. Current Multipoles 

The dominant current multipole moment is 



1(1 + 1] 



5 J, ^ 2lu 



pa 



1+?, 



Kida, 



which is 



a'^lKA^da, 



(23) 



6Ji 



which gives 
(l + lf 



pnR 



1+4: 



(32) 



(33) 



E 



2 n5 



8/ 



pfl^R 



The contribution of he dominant multipole moment to 
(24) dE/dt\g^ is an 0{Q?^+'^) and an 0{Q?^+^) term. 



5.2. Dissipation due to Gravitational Waves 



The dissipation of energy due to the emission of gravita- 
tional waves can be estimated from the standard multipole 
formula: 

dE 



5.2.3. Growth-time 

The growth time due to the emission of gravitational 
waves is 



2E 



dt 

gw 

where 
Ni = 47r 



o,Y.Niu:f+W5DT\^ + \5.n 



(Z + l)(/ + 2) 



m\2\ 



In (25) 



5D'7 



1{1-1)[{21 + IW 



6pa%'^*dV, 



(25) 



(26) 



(27) 



dE/dt\ 



(34) 



gw 



Including both the mass and current multipole contribu- 
tions and keeping the frequency to 0(il'^), we obtain 



c2'+3 ^{l + 1)3 



are the mass multipole moments and 



G il 
([ + 11' 

xMR'^'n^'+^ 



(do + (J2'^^) [l — ao — (72^^^] 



2121 + 1 



Ni+ [I -(To- cTztn^] Ni+im'^ 



5Jr = 2^111 j a' {p5v + 5pv)YT''*dV (28) 

are the current multipole moments, where v is the ve- 
locity of the fluid and Y"^^* are the "magnetic" vector i 
harmonics (see Thorne 1980). 



To lowest order in 51, (p3) reduces to 



[(2? + 3)!!]2 



l + \ 



24G (2; + 3)(;-l)2' \l + 2 



21+2 



-21-2 



MR^^ 



(35) 



(36) 
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5.3. Dissipation due to Shear Viscosity 



6. Critical Angular Velocities 



The dissipation of energy because of the shear viscosity of 
neutron star matter is 
dE 



-2 j rjSa^'^SaltdV, 



dt 
where 



(37) 



(38) 



(see e.g. Ipser &: Lindblom 1991 ) and rj is the shear vis- 
cosity coefficient. We obtain 

a'^\adaKi\^da 



dE 



4{l + l)ujf.T] 



+ il-l){l + 2) / a'\Ki\'da 
Jo 



which yields 
dE 
'dt 



{I + - l){2l + 

4:1 



The damping time due to shear viscosity is 
3 M 



47r(/- l)(2Z + 3) riR' 



(39) 



(40) 



(41) 



5.^. Dissipation due to Bulk Viscosity 

In a neutron star, bulk viscosity can arise, because of 
the departure from nuclear reaction equilibrium, such as 
beta-equilibrium, during the compression and expansion 
of matter caused by an oscillation. The energy is dissi- 
pated at a rate 

'^^ ^-JcM^dV, (42) 



dt 



bv 



where C is the coefficient of bulk viscosity 
Ap 



Sa 



(43) 



is the expansion of the fluid and F is the adiabatic index. 
The last relation follows from baryon conservation in an 
adiabatic oscillation. Strictly speaking, in a uniform den- 
sity star, (5(T = 0. But, we assumed the uniform density 
approximation only to make calculations easier. For the 
bulk viscosity we use an approximate time scale, that has 
been derived by Cutler fc Lindblom 1987 for spheroidal 
oscillations in uniform-density stars. Since the bulk vis- 
cosity arises because of the change in density, for toroidal 
oscillations we use the spheroidal formula, but with I re- 
placed by ^ -I- 1: 

_ 3(2/ H- 5) F^Af 
^''^ ~ 27r(/ + 1)3 CR ■ ^ ^ 

For the purpose of estimating the bulk viscosity only, F 
is taken to be equal to 5, i.e. correspond to that of a stiff 
(nearly uniform density), N = 0.25 polytrope. 



Below the superfluid transition temperature, which is 
T ^ 10^ K, the shear viscosity is dominated by electron- 
etectron scattering and an approximate formula for the 
viscosity coefficient is 



7? = 6x lO^s^ g/cms. 



(45) 



(Cutler & Lindblom 1987), where the notation pis means 
normalization of the density to 10"'^^ gr/cm^ and Tg nor- 
malization of temperature to lO^K. Above the superfluid 
transition temperature, the shear viscosity coefficient due 
to neutron-neutron interactions is 



9/4 

7/ = 2x 10i«^ g/cm s. 



(46) 



(Flowers fc Itoh 1979). The bulk viscosity will be impor- 



tant in hot, newly-born neutron stars, but its coefficient is 
not as certain as the coefficient for shear viscosity. Sawyer 
(1989) estimates the bulk viscosity in neutron star matter, 
assuming that the star cools through the modified URCA 
process and that it is transparent to neutrinos. The coef- 
ficient he obtains is 



C = 6 X lO^^p^gW-'Tg^'g/cms 



(47) 



It has been suggested ( Lai &: Shapiro 1995 ) that for tem- 
peratures larger than a few times Vd^K the neutrino opti- 
cal depth is still large and the bulk viscosity is thus inac- 
tive. If the star cools through the direct URCA reaction, 
the bulk viscosity will be much larger than in (^^, but 
again only for temperatures for which the star is trans- 
parent to neutrinos. It becomes apparent that, depending 
on the cooling process and on the neutrino optical depth 
in a newly-born neutron star, the bulk viscosity can al- 
most completely damp non-axisymmetric instabilities or 
have only a small effect on them. A more detailed study is 
needed and departure from equilibrium of other interac- 
tions (such as interactions between quarks at lower tem- 
peratures) should also be considered. For the time being 
we will use (47) as an conservative average of the large 



error bars associated with the bulk viscosity. 

Another dissipation mechanism that can affect the in- 
stability is the superfluid mutual friction. Estimates by 
MendeU (1991) and Lindblom & MendeU (1995), suggests 
that mutual friction could suppress the gravitational - ra- 
diation - driven instability of /-modes, when the temper- 
ature of the star is between 10^K< T < lO^K. It is not 
clear whether mutual friction will have the same effect for 
r-modes, and a new calculation of this effect is needed. 

At each value of the temperature of the star, the criti- 
cal angular velocity above which gravitations radiation has 
the shortest time-scale, compared to the viscosity time- 
scales, is obtained by solving the equation 



1 



1 1 
— H 

^SV 'HiV 



0. 



(48) 
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We specialize to a specific neutron star model with ra- 
dius R = 12.47km and mass M ~ 1.5Mq (same as 

The 



in Andcrsson ct al. 



1998 and Lindblom et al. 
3.4 X 10"gr/cm^ 



1998) 



density of the star is p ■ 

Fig. 1, shows the critical angular velocity as a func- 
tion of temperature (in units of the angular velocity at 
the mass-shedding limit for Newtonian, uniform density 
and uniformly rotating stars, ^Ik — O-^Ty/irGp). The solid 
curve corresponds to the 0{n'^) equation (|35|). A rapidly- 
rotating neutron star, born at temperatures lO^^K loses 
angular momentum because of the r-mode instability and 
slows-down. The minimum angular velocity it could reach 
is of = 414S-1 (or period of 15ms) at T ~ 1 x 1.5 x 10^ 
K. The mass of the neutron star, or the adiabatic index, 
do not have a significant effect on the minimum critical 
angular velocity. The radius of the neutron star, however, 
(which can range from 10km to 15km), does have a con- 
siderable effect and the radius of our model represents a 
mean value of the expected radius of a typical 1.5Mq neu- 
tron star. 



For low rotation rates, one can use (36) to construct 
approximate equations for the two parts of the curve in 
Fig. I. The part of the critical curve where the shear vis- 
cosity dominates can be approximated by 



ni'^^ = 581 



/I0km\ 



3/2 



V R 



lO^K 
T 



1/3 



(49) 



wihle the bulk viscosity dominated part is described by 



362 



R 



lOkr 



9/8 



/ T 



VlO^K 



3/4 



(50) 



The two approximate expressions are shown as dotted 
curves in Fig. 1. For a period 1.56 ms (the period of the 
fastest known millisecond pulsar), the lowest order critical 
angular velocity differs from the O(il^) result by ~ 17%. 

In a similar way, the lowest order approximations to 
the dissipation timescales are 



i<.w = 1.4 X lO'' 



lO^s 



3„-l 



n 



lAMa\ /lOkmX'^ 



M 



\ R J 



t^^ = 3.6 X 10' 



and 



1AM, 



(lOkm) (iO^k) ^ M 



s, (51) 



(52) 



= 4.6 X 10^ 
/ R 



V 10km 



n 



103 



1.4M, 



M 



lO^K 
T 



(53) 



Our current results for the onset of the r-mode instabil- 
ity correspond to neutron stars with a very stiff equation 
of state. The results in Andersson et al (1998), correspond 
to a much softer equation of state (an N=1.0 polytrope) 
and a comparison is shown in Fig. 2. The minimum crit- 
ical temperature is roughly the same for both equations 



a. 
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Fig. 1. Critical angular velocity for the onset of the r- 
mode instability as a function of temperature (for a 1.5 
Mq neutron star model) . The solid line corresponds to the 
0{fl^) result using superfluid (e~— e~) shear viscosity, and 
Sawyer's (1989) estimate for the bulk viscosity. Dotted 
lines are lowest order approximations, while the dashed 
line corresponds to normal matter (n — n) shear viscosity. 



of state, although it occurs at a somewhat smaller tem- 
perature in the uniform density case. This shows that the 
uncertainty in the equation of state does not have a sig- 
nificant impact on the r-mode instability. 

7. Discussion 

Our analytical results for the onset of the r-mode instabil- 
ity in neutron stars agree well with numerically obtained 
results, for the same neutron star model. Using our ana- 
lytic formulae, the uncertainty in the bulk viscosity can be 
easily explored for different present and future estimates of 
the bulk viscosity coefficient. The shear viscosity for lower 
temperatures is also uncertain, since a high shear viscosity 
due to the mutual friction between superfluid vortices be- 
low the superfluid transition temperature could suppress 
the instability. If future investigations provide a definite 
answer on the effect of mutual friction in a superfiuid, one 
can easily study the implication on the r-mode instability 
using the analytic formulae presented in this paper. 

We would also like to point out that, although an Z = 1 
dipole mode does not radiate in a non-rotating star, it does 
emit gravitational waves through the coupling to higher 
order terms, in rotating stars. According to (plj), the shear 
viscosity for I — 1 vanishes and thus cannot affect the 
emission of gravitational waves, this is not true for real- 
istic equations of state but still the damping times are 
extremelly long. The I — 1 r-mode will radiate through 
the coupling to spheroidal I = 2 terms, i.e. it will generate 
mass quadrupole radiation. The frequency of this mode 
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Fig. 2. Critical angular velocity for the onset of the r- 
mode instability as a function of temperature (for a 1.5 
Mq neutron star model). The present uniform density re- 
sult [N = 0) is compared to the critical curve for the 
equation of state used in Andersson et al. (1998), [N — 1.0 
polytrope). The minimum value of ^d^K is roughly the 
same, and the effect of the equation of state is mainly to 
shift the critical curve to different temperatures. 
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in the rotating frame is ujr = while in the inertial 
frame the frequency is uii = (3/4) nj'^. According to the 
criterion for the onset of the CFS-instability, the I = 1 
r-modes are thus stable to the emission of gravitational 
waves, in contrast to the I > 2 modes. Such stable oscilla- 
tions, unaffected by shear viscosity at low temperatures, 
could be excited during a neutron star glitch. In analogy 
to the / = 1 r-modes, an ^ = 1 spheroidal mode, like the 
/ mode or the p-modes, will emit current quadrupole ra- 
diation and this case needs further study. 

Acknowledgements. We would like to thank N. Andersson, 
J. L. Friedman and B. F. Schutz for helpful discussions and 
suggestions. N.S. acknowledges the generous hospitality of the 
Max-Planck-Institute for Gravitational Physics in Potsdam, 
Germany. 



References 

Andersson, N. 1998 A new class of unstable modes of rotating 



relativistic stars, ApJ in press, gr-qc/9706075 
Andersson, N., Kokkotas, K. D. & Schutz, B. F. 1998, Gravi- 
tational r adiation limit on t he spin of young neutron stars, 
preprint, astro- ph/9805225| 



Chandrasekhar, S. 1970, Phys. Rev. Lett., 24, 611 
Chandrasekhar, S. & Ferrari, V. 1991, Proc. R. Soc. Lond. A, 
433, 423 

Cutler, C. & Lindblom, L. 1987, ApJ, 314, 234 
Flowers, E. and Itoh, N. 1979, ApJ, 230, 847 
Friedman, J. L. & Schutz, B. F., ApJ, 22, 281 



